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Double quantum dot with tunable coupling in an enhancement-mode silicon 
metal-oxide semiconductor device with lateral geometry 
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We present transport measurements of a tunable silicon metal-oxide-semiconductor double quantum dot 
device with lateral geometry. Experimentally extracted gate-to-dot capacitances show that the device is 
largely symmetric under the gate voltages applied. Intriguingly, these gate voltages themselves are not 
symmetric. Comparison with numerical simulations indicates that the applied gate voltages serve to offset 
an intrinsic asymmetry in the physical device. We also show a transition from a large single dot to two well 
isolated coupled dots, where the central gate of the device is used to controllably tune the interdot coupling. 
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Recent progress towards demonstrating electron spin- 
based quantum bits in semiconductor quantum dots 1 3 
has lead to renewed interest in fabrication of quantum 
dot structures in silicon. Silicon is a strong candidate 
due to a relatively small electronic spin-orbit coupling, 
low concentration of nuclear spins (a potential source of 
electron spin decoherence), and the ability to leverage 
mature silicon fabrication technologies. 

We report low-temperature transport measurements of 
a Si metal-oxide semiconductor (MOS) double quantum 
dot (DQD). In contrast to previously reported measure- 
ments of DQD's in Si MOS structures^ - —, our device has a 
lateral gate geometry very similar to that used by Petta 
et al 2 to demonstrate coherent manipulation of single 
electron spins. This gate design^ provides a high degree 
of tunability, allowing for independent control over indi- 
vidual dot occupation and tunnel barriers, as well as the 
ability to use nearby constrictions to sense dot charge 
occupation^. 

Figure 1(a) shows a top-down scanning electron mi- 
croscope (SEM) image of the partially processed device 
and a schematic of the final device cross section. The 
upper Al top gate is used to accumulate carriers at the 
Si-Si02 interface. The lower, patterned polycrystalline 
silicon (polysilicon) gates are used to deplete carriers to 
define the DQD region. The device is fabricated on a 
lightly p-type Si substrate (2 - 20 Q cm). Initial fabri- 
cation steps (Ohmic contact formation to channel, gate 
oxide growth, and polysilicon deposition and doping) are 
performed in a fully qualified CMOS facility. Ohmic con- 
tacts (not shown in Fig. 1(a)) are formed by implantation 
of As, followed by a 900 °C, 15 min. activation anneal. 
Next, a 35 nm Si02 gate oxide is grown in dry O2 at 
900° C, with a subsequent N2 anneal at 900° C for 30 min., 
immediately followed by polysilicon deposition (200 nm) 
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and doping. The lower-level polysilicon depletion gates 
are defined by e-beam lithography and dry etch and are 
200 nm tall with a ~100 nm linewidth for the narrowest 
features (see Fig. 1(a)). After the polysilicon etch, the 
device is exposed to a thermal oxidation step which forms 
a 30 nm layer of Si02 on the surface of the polysilicon 
gates. Finally, the polysilicon is further isolated from the 
global Al top gate by 60 nm of AI2O3 formed via atomic 
layer deposition (ALD). More details about the fabrica- 
tion of this device and its operation in single-dot mode 
can be found in Ref. [l2|. 

The device conductance is experimentally determined 
via standard low-frequency lock-in measurements with an 
rms ac source-drain bias of 10 - 50 uV. Unless otherwise 
noted, the Al top gate TG, T, L, and R (see Fig 1(a)) 
voltages are held constant at Vtg = 5 V, Vt = -0.3 V, 
Vl = V, and Vr = -2 V. All measurements shown are 
performed in a dilution refrigerator with a temperature 
of T ~ 20 mK. We note that although the conductance of 
the device continues to evolve below a fridge temperature 
of 100 mK, we cannot be certain our dot electron temper- 
ature is equal to our fridge temperature in this regime. 
However, the precise value of the electron temperature 
should not affect the conclusions of this letter. 

Figures 1(b)- (e) show DQD conductance versus left 
and right plunger gate voltages Vlp and Vrp at four dif- 
ferent center plunger voltages Vcp, demonstrating the 
ability to tune the DQD from a highly-coupled regime, 
where the transport is reminiscent of that expected for a 
large, single dot, to a weakly-coupled regime, where con- 
duction can only take place at the so-called triple points. 
As expected, the interdot coupling shows a strong de- 
pendence on Vqp, whereas we find a nearly negligible 
dependence on Vlp and Vrp. 

In Fig. 2(a) we show DQD conductance versus Vlp 
and Vrp at Vcp = -1-2 V. The capacitance between the 
plunger gates and dots can be extracted from the dimen- 
sions shown in Fig. 2(a): Clp-ld = e/5VLp_LD = 8.1 




FIG. 1. (a) SEM image of partially processed device showing polysilicon gates and schematic of final device cross-section, (b) - 
(e) Conductance of double dot G versus Vlp and Vrp, showing transition from a large single-dot to a well-defined double dot, 
at CP gate voltages (b) V C p = -0.6 V, (c) V C p = -0.8 V, (d) Vcp = -1.0 V, (e) V C p = -1.2 V, for V T = -0.3 V, V L = V, and 
V R = -2 V. 
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FIG. 2. (a) DQD conductance G vs. Vlp and Vrp at Vcp = 
-1.2 V, showing LP (RP) voltage required to change the occu- 
pation of the left dot (right dot) by 5 (4) electrons: 5£Vlp-ld 
(4£Vrp_rd)- (b) Closer view of triple points at Vcp = -1.2 V 
for zero dc source drain bias, (c) Triple points at Vcp = -1.2 V 
for dc source-drain bias V s d = 0.5 mV, showing dimensions of 
bias triangles, (d) Interdot coupling F versus Vcp determined 
from plunger gate voltage separation of triple points. 



aF, Crp_rd = e/^^RP-RD = 5.3 aF. Figure 2(c) shows 
bias triangles at a finite dc source-drain bias of 0.5 mV 
at Vcp = -1.2 V. Using C lp _ ld( rp_r D) /C ld(rd) = 
V r s d/AVi J p(Rp) we find total dot capacitances Cld = 210 
aF and C RD = 170 aE±i 

Figure 2(b) shows a closer view of single pair of triple 
points at Vcp = -1.2 V and zero dc source-drain bias. 
The distance between triple points in plunger gate volt- 



age space SVs, as shown in Fig. 2(b), can be used to 
estimate the interdot coupling. In Fig. 2(d) we show the 
fractional splitting of the triple point F versus CP gate 
voltage, where F = 2SVs/SVp (see Fig. 2(a), (b)) is a 
measure of the interdot coupling, defined as the ratio of 
the diagonal separation between triple points to the sep- 
aration between charge domain a 14 i 15 . The error bars are 
determined by the ability to visually resolve the position 
of the triple points on a honeycomb plot. Figure 2(d) 
shows that sweeping CP smoothly varies F from nearly 
zero to one, demonstrating the ability to tune the DQD 
from weakly-coupled to a regime where the dots are fully 
merge d 14 ! 16 . 

Figure 3(a) shows a comparison between experimen- 
tally determined capacitances from the various gates to 
the dots and values determined via modeling. The val- 
ues Cmeas in Fig. 3(a) are obtained from honeycomb 
plot dimensions from double dot transport, or from the 
gate voltage period of Coulomb blockade oscillations in 
single dot transport. For double dot transport, the gate 
voltages are Vtg = 5 V, V T = -0.3 V, V L = V, Vi P 
= V, Vcp = -1-2 V, Vr P = V, V R = -2 V. Values 
marked with an asterisk indicate capacitances from sin- 
gle dot transport, obtained for the left dot by setting Vr 
= +1 V, or for the right dot by using Vl = +1 V, while 
maintaining all other gate voltages the same as used for 
double dot transport. Error bars are determined by vari- 
ations in the charge transition period. 

Figure 3(b) shows a contour plot of calculated elec- 
tron density in the DQD region at gate voltages equal 
to those used for the experimental capacitances. The 
calculations are semiclassical (Thomas- Fermi approxima- 
tion) and were performed using TCAD Sentaurusii, a 
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FIG. 3. (a) Table of capacitances between dots (right dot 
(RD) and left dot (LD)) and gates (as labeled in Fig. 1(a), 
comparing measured and calculated values. Cmeas are ex- 
perimental values for the case Vcp = -1.2 V, Csentaurus val- 
ues are computed using a semiclassical model for dot electron 
density, and Ccfd-ace values are calculated by treating the 
dots as metallic regions defined by the shaded areas sketched 
in (b). Asterisks indicate capacitances from single-dot trans- 
port, obtained by setting Vr (Vl) = +1 V for left (right) dot 
measurement, (b) Calculated electron density in dot regions 
relative to poly silicon gates, where the dotted lines are con- 
tours of electron density from 1 to 7 x 10 11 cm -2 in 1 x 10 11 
cm -2 increments. The dashed lines define the boundary be- 
tween the left and right dot regions and between the dot re- 
gions and the leads, (c) Sketch of dot region used to obtain 
good agreement between experimental and calculated capac- 
itances. Shaded area represents dot region, which is treated 
as a metallic (perfect conductor) region in order to calculate 
capacitances. 



not reflect transport in the actual device, which yields 
C*r-rd = 11-6 aF versus Cl-ld = 17.3 aF. The dot re- 
gions shown in Fig. 3(c), which result in reasonable agree- 
ment between Cmeas and Ccfd-ace, predict a right dot 
that is smaller than the left dot, but not as dramatically 
so as suggested by the contours in Fig. 3(b). 

The disagreement between Cmeas and Csentaums sug- 
gests that the Sentaurus calculation fails to capture some 
aspect of the actual device, such as unintentional imper- 
fections that may be present in the experiment, but are 
not included in the calculation. For example, a spatially 
varying distribution of charge in the dielectric regions or 
variations in film thickness (e.g. non-uniform AI2O3 de- 
position) can lead to unintentional variations in electron 
density in the dot region. However, the relative symme- 
try of the dot contours shown in Fig. 3(c) (as compared 
to Fig. 3(b)) suggests that by use of asymmetric gate 
voltages (Vr = — 2 V versus Vl = V) we are able to 
partially compensate for this disorder. 

In conclusion, we have demonstrated a gate-defined 
Si MOS double quantum dot with lateral geometry. 
Comparison between experimental and calculated capac- 
itances suggest the presence of disorder and that we are 
able to partially compensate for this disorder by adjust- 
ment of gate voltages, due to the tunability of our lateral 
geometry. Our data also shows the ability to control the 
interdot coupling over a wide range of energies. Similar 
device structures may allow for further study of DQDs 
in the few electron regime, which would be of interest for 
quantum computing applications. 



commercial device simulation package. The capacitances 
Csentaurus are obtained by integrating the charge density 
over the entire left or right dot region and calculating 
the change in the integrated dot charge due to a small 
change in gate voltage. The boundaries used to separate 
the left dot from the right dot and the dot region from 
the leads are shown by the dashed lines in Fig. 3(b). 

In Fig. 3(c), we show a model dot region (shaded 
area) used to calculate capacitances between the various 
gates and dots by treating the dot region as a perfectly- 
conducting metallic sheet. These capacitance calcu- 
lations were performed using CFD-ACE+ 18 , a finite- 
element modeling package. The shape of the dot was 
chosen by starting with a region which approximately 
follows the contours defined by the gates, and then mod- 
ifying the distance between the dot and various gates as 
to obtain good agreement with the experimental capaci- 
tances Cmeas listed in Fig. 3(a). 

The agreement between the capacitances predicted by 
the Sentaurus calculation and those obtained from ex- 
periment is relatively poor. As shown in Fig. 3(b), the 
Sentaurus calculation predicts an asymmetric dot region, 
resulting in a capacitance of Cr_rd = 4.6 aF from 
the R gate to right dot, versus Cl-ld = 34.4 aF from 
the L gate to left dot. This extreme asymmetry does 
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